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ABSTRACT 

The Local Group is a small galaxy cluster with the membership of 62 nearby galaxies 
including the Milky Way and M31. Although the Local Group has yet to be virial- 
ized, it interacts with the surrounding matter as one gravitationally bound system. 
To understand the formation and evolution of the Local Group as well as its member 
galaxies, it is important to reconstruct the gravitational potential field from the sur- 
rounding matter distribution in the local cosmic web. By measuring the anisotropy in 
the spatial distribution of the Local Group galaxies, which is assumed to be induced 
by the local gravitational tidal field , we resolve the three dimensional structure of the 
gravitational potential in the vicinity of the Milky Way smoothed on the Local Group 
mass scale. Our results show that (i) the minor principal axis of the Local Group 
tidal field is in the equatorial direction of a p = 15 h 00 m and S p = 20 d ; (ii) it has a 
prolate shape with axial ratio of 0.5 ± 0.13; (iii) the global tides in the Local Group is 
quite strong, which may provide a partial explanation for the low abundance of dwarf 
galaxies in the Local Group. 

Key words: methods:statistical - cosmology:theory - galaxiesxlustering - galax- 
ies:halos - large-scale structure of Universe 



1 INTRODUCTION 

Since iHubbld 1 19361 ) first recognized its existence, the Lo- 
cal Group has been wid ely investigated by various observa- 
tions (for a review, see lMateolll99Sl ; IVan den Berghl l200l 
and references therein). It is now well known that the Lo- 
cal Group consists of two massive halos (the Milky Way 
and M31), their satellites, and the neighboring galaxies. It 
has been estimated by numerical simulations based on the 
disk formation model that the dark matter halos of the 
Milky Way and M31 have virial radiu s of approximately 
200-300 kpc separated by a bout 700 kpc (jKlypin et alj|2002i; 
IWidrow fc Dubinski|[2OT5l ). The comparison of the numeri- 
cal results with observational constraints suggested that the 
Local Group is not yet relaxed with two main halos in the 
merging process. 

Despite that the Local Group has yet to be virialized, 
it is also true that the Local Group behaves as one gravita- 
tionally bound system when i nteracting with the s urround- 
ing large-scale structure (e.g. JVan den Bcrgh 200fJ). There- 
fore, to understand the formation and evolution of the Local 
Group under the gravitational influence from the surround- 
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ing matter distribution, it should be desirable to resolve the 
gravitational potential field around the Local Group. 

As pointed out bv lSpringel et alj l|2004j ) , it is the struc- 
ture of the gravitational potential that reflects more directly 
the three dimens ional distribution of d ark matter in a given 
region. Recently. Irlavashi et al.l (120071 ) showed by numerical 
analysis of the N-body simulation results that the poten- 
tial field is much smoother than the density field and the 
iso-potential surfaces are well approximated by concentric 
triaxial ellipsoids. 

We attempt here to resolve the three dimensional struc- 
tural shape of the gravitational potential in the regions 
around the Local Group, assuming that its iso-potential sur- 
faces can be well approximated as a triaxial ellipsoid. To 
achieve this goal, we emplo y the analytic algorithm devel- 
oped bv lLee fc Kand (|2006l . hereafter, LK06) which allows 
us to determine the eigenvalues of the tidal tensor (defined 
as the second derivative of the gravitational potential) from 
the anisotropic spatial distribution of the galaxies in a given 
cluster. 

The underlying logic of the LK06 algorithm is as follows: 
The cluster galaxies are observed to be preferen tially located 
near the major axes of the parent clusters (|Knebe et al.l 
|2004| ; [Zentner et al.ll2005l ). Recent numerical analyses based 
on N-body simulations have shown that this anisotropic dis- 
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tribution of the cluster galaxies originate from th e local 
tidal field in the p arent clusters (|Atlav et alj|2006h . What 
iLee fc Kand (|2006h found is that it is possible to determine 
the eigenvalues of the local tidal field by measuring the 
anisotropy in the spatial distribution of the cluster galax- 
ies. Once the eigenvalues of the tidal field are determined, 
the axial r atios of the iso-potential surfaces can be readily 
calculated (|Bardeen et al.lll986l ; llond fc Mverslll996r i. 

The plan of this paper is as follows. In §2, the LK06 
model is briefly overviewed and our theoretical formula 
based on the LK06 model is presented. In §3 the observa- 
tional data are described and the anisotropy in the spatial 
distribution of the Local Group member galaxies is mea- 
sured. Also the structural shape of the gravitational poten- 
tial in the region around the Local Group is resolved. In §4 
the results are discussed and a final conclusion is drawn. 



2 THEORETICAL BACKGROUND 

2.1 The gravitational potential profile 

Let < I > (r) represent the gravitational potential field smoothed 
on the Local Group mass scale. In the neighborhood of the 
Local Group, the Taylor expansion gives 

$(r) » $(0) + r^lo + -nr^o^lo, (1) 

where it is assumed that the center of the mass of the Local 
Group is positioned at the origin of r = 0. Supposing that 
the Local Group represents a local peak of the smoothed po- 
tential field, the second term proportional to di&\o in equa- 
tion fl]) becomes zero. The tidal shear tensor T in the Local 
Group is defined as Ty = didj<&\o at the local peak to which 
the third term in equation (fTJ is proportional. 

In the frame of the principal axes of the tidal shear 
tensor, equation |T} can be expressed as 

*(r) - *(0) « ir?Ai, (2) 

where, Ai, A2, A3 are the eigenvalues of T with Ai ^ A2 ^ 
A3. Let ai, ei2, ei3 be the three semi-axes of an ellipsoidal 
iso-potential surface of ^(r) = const. = (f>. By equation @, 
one can relate {eii}? to as 



2(*(0) - <P) 



A, 



1/2 



(3) 




Therefore, the two axial ratios of the ellipsoidal iso-potential 
surfaces can be written only in terms of {Ai}f =1 as 



(4) 



Equation Q implies that once the eigenvalues of the tidal 
shear tensor at the Local Group region are determined, the 
three dimensional shape of the Local Group potential can 
be reconstructed. 



2.2 Overview of the LK06 algorithm 

Assuming that the tidal field in the host halo is responsible 
for the aniso t ropic spatial distribution of the substructures, 
I Lee fc Kand ( 2006) have derived an analytic expression for 



the probability distribution of the cosines of the polar an- 
gles of the substructure position vectors in the tidal shear 
principal axis frame as 



p(cos0) = inti( 1 - s + 3sA I 2 ) 



-1/2 



cos <j> 1 sin sin (f> . cos 9 
s+3sA? l-s+lisX?, l-s+3sA; 



-3/2 



d<t>, (5) 



where s is a correlation parameter in the range of [—1, 1] 
and {Ai}f =1 are the unit eigenvalues of T satisfying the 
constraint of A 2 + A| + A3 = 1. Here, the free parameter 
s quantifies the degree of the alignment between the minor 
principal axis of the tidal tensor and the position vectors 
of the substructures: If s = — 1, then the substructures are 
most strongly aligned with the minor principal axis of the 
tidal tensor; If s — 1, they are most strongly anti-aligned 
while the case of s — corresponds to no alignment. 

Since the minor principal axis of the tidal tensor at the 
region of a given host halo corresponds to the direction in 
which the tidal fo rce is least strong a nd thus the host halo 
is most elongated, ILee fc Kand (2006) assumed that the mi- 
nor principal axis of the tidal tensor is parallel to the major 
principal axis of the inertia tensor of the host halo. Under 
this assumption, they suggested that by measuring p(cos#) 
from the spatial distribution of the substructure positions 
and fitting it to equation ©, one can determine the unit 
eigenvalues of the tidal tensor and the value of the corre- 
lation parameter, which will in turn yield the values of the 
axial ratios of t he host halo. To de monstrate the validity of 
their algorithm. ILee fc Kand l|2006h tested it against N-body 
results, and showed that their algorithm indeed reconstruct 
the actual shapes of cluster halos within 20% error. 

In the original approach, ILee fc Kand (2006) used the 
virialization constraint of Ai = 8 C wher e S c ~ 1.68 is th e 
density threshold for a virialized structure l|Eke et al.lll996T l . 
However, note that the LK06 algorithm can be applied even 
to a region that is not completely virialized, since the virial- 
ization condition was not used to derived equation ([5]) itself. 
In fact, it can be applied to any region as long as the member 
objects in the region show tidally induced anisotropy. For 
instance, it can be used to describe the anisotropic spatial 
distribution of galaxy clusters in the superclusters due to the 
supercluster tidal field (i.e., the linear tidal field smoothed 
on the supercluster mass scale) , even though the su perclus- 
ters are not virialized objects (ILee fc Evrardl [20071 ). It has 
been also used to describe the anisotropic spatial distribu- 
tion of void galaxies due to the tidal effect from the sur- 
rounding f ilaments even tho ugh the voids are obviously not 
virialized (ILee fc ParkH2006T ). 

When the LK06 algorithm is applied to an unviriazlied 
structure, the additional constraint of ^ 3 Ai = 8 C should 
not be used. Therefore, what one can reconstruct by the 
LK06 algorithm for this case is not the shape defined by 
the inertia momentum tensor of the structure but the shape 
of the potential around the structure. In the following sec- 
tion, we reconstruct the three dimensional structural shape 
of the gravitational potential in the region around the Local 
Group by measuring the values of s and {Ai}f = i from recent 
observational data. 
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Figure 1. Probability density distribution of the cosines of the 
angles between the position vectors of the Local Group member 
galaxies in the CL frame and four different choices of major prin- 
cipal axis of the Local Group. In each panel, the histogram with 
Poisson errors represents the observational data points of Local 
Group, the solid line is the analytic fitting from LK06 recon- 
structing algorithm, and the dotted line corresponds to the case 
of no correlation. Four plots in each panel show how the agree- 
ments between the observational and the analytic results change 
if different directions are used rather than the major principal 
axis. 



3 APPLICATION TO THE LOCAL GROUP 
DATA 

3.1 Tidal Alignment of the Local Group 

We use the observational data of 62 member galaxies be- 
longed to the Local Group. Table [1] lists the latest infor- 
mation on the equatorial coordinate (a, 5), distance (r) in 
unit of Mpc and absolute magnitude (M v ) of each member 
galaxy. The information on a, S, r and M v for the #1 — 40 
and # 48 galaxies are obtained from iKarachentsev et al.l 
|200J), and for the other galaxies, we show in the seventh 
columns the references from which we obtained the informa- 
tion. For those galaxies marked by ★ in Table [T] we used th e 
updated information given in iMcConnachie et al.l l|200ot) . 
Since the mass information is not available, we calculate 
the center-of-luminosity (CL) instead of the center-of-mass 
(CM) for the Local Group by using the information on the 
absolute magnitude of each member galaxies. The Cartesian 
position of the Local Group CL in the equatorial coordinate 
system is found to be CL = [0.351, 0.069, 0.259] in unit of 
Mpc. 

We measure the three dimensional positions of the Lo- 
cal Group member galaxies in the CL frame in the equa- 
torial coordinates. Then, we attempt to find a direction in 
the equatorial coordinates, say, a p and S p , with which the 
position vectors of the Local Group member galaxies in the 
CL frame show the highest degree of anisotropy. We will re- 



X 

cn 
o 




0.5 - 



100 200 300 
a. [deg] 




-50 



50 



<5 P [deg] 



Figure 2. The variation of x as a 



and 



5 P (bottom) when the other parameters are fixed at the best-fit 
values. In each panel the vertical dotted line indicates the position 
at which x 2 has the minimum value (a v = 225°, S p = 20°). 



gard this direction, if found, as the direction of the minor 
principal axis of the local tidal field. 

Let 6 represent the angle between the equatorial po- 
sition vector of a given member galaxy and the direction 
of (a p , S p ). Changing the values of a p and 8 P by 5° sys- 
tematically from 0° to 360° and —90° to 90°, respectively, 
we measure repeatedly the probability distribution, p(cos6>), 
and search for (a p , S p ) which leads to the maximum devia- 
tion of p(cos(9) from uniform distribution. The bin size of 5° 
is chosen as the smallest bin size below which the numerical 
noise dominate. By fitting equation ((5]l to the observational 
results, we determine a direction of (a p , S p ) which leads the 
maximum deviation of p(cos 6) from uniform distribution. 

Fig. [T] plots the cases of four different directions of the 
minor principal axis of the Local Group tidal field as exam- 
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Table 1. Observational Data on the Local Group 
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Figure 3. The mean values, (a p ) (top) and (S p ) (bottom), aver- 
aged over 10 samples of N randomly selected galaxies out of the 
62 Local Group members. In each panel the errors represent one 
standard deviation between samples and the horizontal dotted 
line correspond to our result of (a p = 225°, S p = 20°). 



pies. In each panel, the histogram with Poisson errors cor- 
responds to the observational results while the solid curve 
is the fitting model (eq.[5]) with the best-fit values of Ai, 
A2 and s: For a given direction of (a p , 5 P ), we fit equation 
((5| to the observational result of p(cos#) by adjusting the 
values of Ai, A2, and s and determine for the best-fit values 
of the three parameters with the help of the x 2 statistics. As 
can be seen in Fig. [T] the observed distribution of p(cos9) 
behaves quite differently as the direction of (a p , S p ), changes. 
For a certain direction of (a p ,S p ), the distribution becomes 
quite uniform (top-left). While for another direction p(cos 9) 



increases sharply with (cos 9) indicating preferential align- 
ments of the galaxy's position vectors with the given direc- 
tion and the analytic fits work well (bottom-right). 

It is found that the direction of (a p = 15 h 00 m , 5 P = 
20 d ) leads to the maximum alignments. In other words, when 
this direction is chosen as the minor principal axis of the 
Local Group tidal field, equation © agrees with the obser- 
vational results best with the value of x 2 minimized. Fig. 
[2] plots x 2 as a function of a p (top) and 5 P (bottom) when 
the three other parameters are fixed at their best-fit values. 
In each panel the vertical dotted line indicates the position 
at which x 2 has the minimum value. It is worth mentioning 
here that the direction we found as the minor principal axis 
of the Local Group tidal field (a = 15 h 00 m , S = 20 d ) does 
not coincide with the direction parallel to the connection 
line between M31 and MW. It indicates that the observed 
anisotropic distribution of the LG galaxies is induced by the 
tidal effect of the surrounding matter distribution not by 
the fact that the majority of them are the satellites of M31 
and MW. 

To examine how stable our result is, we construct 10 
samples of N randomly selected galaxies out of the parent 
sample of 62 members. For each sample, we determine the 
best-fit values of a p and 8 P using the same procedure de- 
scribed in the above, and calculate the means of a p and 5 P 
averaged over 10 samples. Fig. 3 plots (a p ) and (6 P ) as a 
function of N as histogram in the top and bottom panel, 
respectively. In each panel the errors represent one stan- 
dard deviation between samples and the horizontal dotted 
line represents our best-fit result (a p = 225°, S p — 20°), 
obtained from the 62 LG galaxies. As can be seen, the vari- 
ations of (otp) and (<5 P ) with N stay within one standard 
deviation, until the value of N decreases down to 45. At the 
bin of N = 40, however, (8 P ) suddenly drops by more than 
two standard deviations. Therefore, one can say that our re- 
sult is stable as long as more than 45 LG galaxies are used 
to reconstruct the minor principal axis of the tidal field. 

Since the effect of the nonlinear process like merging 
between MW and M31 should be dominant near the LG 
center, the spatial distribution of those galaxies located near 
the LG center like the satellites of the MW and M31 are 
not good indicators of the tidal effect. Only when a large 
fraction of distant member galaxies whose positions are not 
severely modified by the merging effect are included, the 
overall spatial distribution of the LG member galaxies will 
reveal the anisotropy induced by the local tidal effect. 

When the minor principal axis of the Local Group tidal 
field is determined as this direction, the values of the fitting 
parameters are found to be Ai = 0.7 + 0.04, A 2 = 0.69 + 0.04 
and s = —0.99 + 0.27. Note that this corresponds to the unit 
eigenvalues of the tidal field smoothed on the Local Group 
mass scale. Fig. [5] compares this best-fit model (solid line) 
with the observational result (dots) with Poisson errors for 
the minor axis direction of (a p = l^OO" 1 , 5 P = 20 d ). As 
can be seen, the fitting model agrees with the observational 
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Figure 4. 65%, 95% and 99% contour plots for s and Ai (left) and for s and A2 (right). In each panel the other parameter is fixed at 
its best-fit value. 



result pretty well. The probability distribution of p(cosO) 
with the best-fit parameters indeed shows strong alignments 
between the minor axis of the Local Group tidal field and 
the positions of the Local Group galaxies. 

The unmarginalized error in the measurements of the 
best-fit values of Ai, A2 and s are calculated b y using the for- 
mula given in l|Bevington fc Robinson! [l 996). For instance, 
the error in the measurement of Ai is calculated as 



AA lv /2(x!-2xi+x|)- 



(6) 



where AAi is the bin size of Ai. Here, \i ~ X 2 (^i»)|i=i,2,3 
and A12 = An + AAi, A13 = A12+AA1 and An is the best-fit 
value. The error associated with the best-fit values of Ai and 
s are also calculated in a similar manner. To account for the 
parameter degeneracy and estimate the marginalized errors 
in the measurement of the best-fit values, we also compute 
the 68%, 95% and 99% confidence regions for (s, Ai) and 
for (s, A2), which are shown in in the left and right panels 
of Fig. |4j respectively. 

By expressing equation Q in terms of the unit eigen- 
values, we can now measures the three dimensional shape of 
the gravitational potential in the region around the Local 
Group: 



ai/03 



02/03 



(7) 



The axial ratios of the ellipsoidal iso-potential surfaces are 
calculated as 0,1/0,3 = 0.51 ± 0.13 and 02/03 = 0.52 ± 0.13. 
The associated errors are estimated by means of the error 
propagation method l)Bevington fc Robinsonlll996h as 



2 2 



df_ 

dXi 



■2 



df_ 

d\ 2 



(8) 
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Figure 5. Probability distribution of the cosines of the angles 
between the spatial positions of the Local Group member galax- 
ies and the major principal axis, a = 15 h 00™ , 6 = 20 d . The solid 
line corresponds to the fitting model (eq. [5]) with the best-fit pa- 
rameters of s = —0.99 and ai/a.3 = 0.51, 02/13 = 0.52 while the 
solid dots with Poisson errors represent the observational results. 



where / represents the two axial ratios. Our results indicate 



that the gravitational potential in the Local Group has a 
prolate shape. Figure [5] plots our final result. 
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3.2 Strength of the Tidal Effect in the Local 
Group 

Now that we have determined the unit eigenvalues of the 
tidal shear field in the region around the Local Group, we 
can estimate the strength of the global tides in the Local 
Group. Since the global tides plays a role of stripping or 
destroying the substructures, if it is found to be strong in 
the Local Group, then it might explain at least partially the 
low abundance of member galaxies in the Local Group. 

The degree of the tidal effect can be quantified in terms 
of the diff erences between the three eigenvalues of the tidal 
tensor as (|Chiueh fc Leell200ll ) 

At = J i [(Ai - A 2 ) 2 + (A 2 - A3) 2 + (A-i - A3) 2 ] . (9) 

Here the value of At is in the range of [0, 1]. The case of 
At = corresponds to no tidal effect. The higher the degree 
of the deviation of At from zero is, the stronger the tidal 
effect is. It is found that At = 0.42 for the Local Group 
from the best-fit values of Ai and A 2 obtained in §3.1, which 
suggests that the global tides in the Local Group is quite 
strong. 



4 DISCUSSION AND CONCLUSION 

We have shown here that the anisotropic spatial distribution 
of the Local Group galaxies can be used to reconstruct by 
means of the LK06 algorithm the shape of the gravitational 
potential in the vicinity of the MW smoothed on the Lo- 
cal Group mass scale and measure the strength of the tidal 
effect, which is important to understand the interaction of 
the Local Group as a whole with the surrounding matter 
distribution and its evolution. 

In their original approach, iLee fc Kand (|2006l ) recon- 
structed the ellipsoidal shape of the host halos rather than 
the gravitational potential from the anisotropic spatial dis- 
tribution of the halo substructu res by using the Zel'dovich 
approximation (Ze rdovichlll97Gi ) as well as the virialization 
condition. We did not apply their original approach to the 
Local Gro up since the Local Group has yet to be complet ely 
virialized (jKlypin et al.ll2002l ; Iwidrow fc Dubinskill2005f ). 

Nevertheless, the LK06 algorithm can in fact apply to 
any overdense region in which the embedded structures show 
anisotropic spatial distribution since it has no t been derived 
under a virialization condition. For example, ILee fc Evrardl 
l|2007l ) used the LK06 algorithm to reconstruct the filamen- 
tary shapes of the superclusters by measuring the cluster- 
supercluster alignments even though the superclusters are 
not completely relaxed systems. For the case of unrelaxed 
systems, however, what can be reconstructed by the LK06 
algorithm is the anisotropy in the gravitational potential but 
not the shape of a bound halo. 

Here, we have reconstructed the anisotropy in the grav- 
itational potential around the Local Group from the spatial 
distribution of the Local Group galaxies. It is shown that the 
Local Group potential is prolate with axial ratios of 0.5±0.13 
with the major principal axis parallel to the direction of the 
equatorial coordinate, a = 15 h 00 m , 6 = 20 d . Even though 
our result still suffers from large errors due to small number 
statistic, it has been obtained without resorting to any ad- 



ditional assumption on the dynamical state or the shape of 
the Local Group. 

It is worth discussing one caveat that our result is sub- 
ject to. As mentioned in Section 3.1, the reconstruction pro- 
cedure is stable as long as more than 40 galaxies of the Local 
Group are used. We think that this caveat comes from the 
radial dependence of the final result. Since the effect of such 
nonlinear process as gravitational merging and interaction 
between the MW and the M31 should be dominant in the 
inner region of the Local Group, the spatial distribution of 
those nearby galaxies like the satellites of the MW and M31 
should not be good indicators of the tidal effect. 

Meanwhile the positions of the distant galaxies located 
in the outer region of the Local Group are not severely mod- 
ified by the nonlinear effect and thus their spatial distribu- 
tion will reflect better the effect of the Local Group tidal 
field. It would be desirable to investigate how the triaxial 
structure of the Local Group potential changes with radial 
distance. It is, however, beyond the scope of this paper since 
the total number of the Local Group member galaxies is too 
small to obtain reliable statistical result on the radial depen- 
dence. We plan to apply the LK06 algorithm to the Virgo 
cluster and investigate the radial dependence of the Virgo 
cluster potential shape, given the fact that the Virgo cluster 
has much more member galaxies than the Local Group. We 
wish to report the result in the near future. 

From the measured eigenvalues of the tidal tensor, we 
have also shown that the the global tides are quite strong 
in the current Local Group system, which may provide an 

explanation to the low abundance of dwarf galaxies in the 

■ 1 1 — 1 

Local Group Ijsomerville! 120021 ) : In the Local Group system 
where the two main galaxies are in the process of filamentary 
merging the low-mass dwarf galaxies are hard to survive due 
to the strong global tides. An interesting question to ask is 
what the critical mass of the member galaxy is to survive 
the global tides in the Local Group? We also intend to work 
on this issue in the future. 
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